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Abstract 

A new method for the measurement of the muon flux in the deep-sea ANTARES 
neutrino telescope and its dependence on the depth is presented. The method is 
based on the observation of coincidence signals in adjacent storeys of the detector. 
This yields an energy threshold of about 4 GeV. The main sources of optical back- 
ground arc the decay of ^''K and the bioluminescence in the sea water. The 
background is used to calibrate the efficiency of the photo-multiplier tubes. 
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1. Introduction 



The ANTARES neutrino telescope is installed in the Mediterranean Sea, 40 km 
offshore from Toulon (France) at a depth of 2475 m [H, [2|. The elements of the 
detector are arranged in twelve vertical lines, each of which is attached to the seabed 
with an anchor and kept taut by a buoy. A detector line comprises 25 storeys 
distributed along the length of the 450 m long electro-mechanical cable, starting 
100 m above the seabed. The storeys are separated by 14.5 m and contain a triplet of 
large-area hemispherical 10" photo-multiplier tubes (PMTs), arranged as in Figure [1] 
0, 0, Isj. The detector has been operated in partial configurations since March 2006 
and was completed in May 2008 j2| . First results are published in [g], Q, @] . 

A large flux of muons produced in the atmosphere by high energy cosmic rays 
passes through the detector. The Cherenkov light produced in water by a muon with 
an energy in excess of 4 GeV can generate correlated signals in adjacent storeys. 
Background light, mainly due to decays and bioluminescent organisms, is also 
present ^]. Although a single decay will produce a relatively small number 
of Cherenkov photons, it can be observed as a coincidence between neighbouring 
PMTs within a storey if the decay occurs in the vicinity. Bioluminescence is a single 



photon process and contributes only to the accidental coincidence rate [10|, lU . 

The present analysis employs a coincidence method to identify decays and 
low-energy atmospheric muons on a statistical basis. The flux of atmospheric muons 
is measured at 24 different depths. The signal is used to calibrate the efficiency 
of the PMTs. 

The paper is organized as follows. The calibration procedure is introduced in 
Section [21 The method to determine the muon flux is presented in Section [3l The 
depth dependence is discussed in Section HJ 



2. Calibration with potassium-40 

Potassium-40 is a radioactive isotope naturally present in the sea water. The 
decay — e~ Vf, ^°Ca yields an electron with an energy up to 1.3 MeV. This 
energy exceeds the Cherenkov threshold for electrons in sea water (0.25 MeV), and 
is sufficient to produce up to 150 Cherenkov photons. Another source of electrons 
above the Cherenkov threshold is the Compton scattering of 1.46 MeV photons, 
which are produced in the process of electron capture -|- e~ — ^°Ar* -|- Ue, 
followed by "^^Ar* ^ "^^Ar + 7. If the decay occurs within a few meters of a detector 
storey, a coincident signal may be recorded by two of the three PMTs on the storey. 
This is referred to as a local coincidence. An example of the distribution of the 
measured time differences between hits in neighbouring PMTs in the same storey is 
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Figure 1: Schematic view of one of the detector storeys. The storey contains three transparent 
glass vessels (17" diameter), each housing a 10" PMT. The PMTs are oriented downwards at an 
angle of 45° from the vertical and with an angle of 120° from each other. The distance between 
the centers of the photo-cathodes of the PMTs is 1.0 m. The storey also contains an electronics 
module for the readout of the PMTs. 



shown in Figure [2] (left). A clear peak around ns is visible. This peak is mainly 
due to single decays producing coincident signals. The data has been fitted 
to a sum of a Gaussian distribution and a flat background. The full width half 
maximum of the Gaussian function is about 9 ns. This width is mainly due to 
the spatial distribution of the '*°K decays around the storey. The rate of genuine 
coincidences is given by the integral under the peak and can be extracted from 
the distribution by subtracting the fiat background of random coincidences. The 
positions of the peaks of the time distributions for different pairs of PMTs in the 
same storey are consistent with zero and are used to cross-check the time offsets 
computed by the time calibration system • 

The spatial uniformity and temporal stability of the rate of decays is used 
to calibrate the relative efficiencies of the PMTs. For every detector storey three 
coincidence rates, r^, and rsi, are measured, which are related to the efficiencies 
of the three PMTs, si, S2 and S3, where the subscript refers to the position of 
the PMT inside the storey. The efficiencies are derived from the measured rates 
assuming that the rate is proportional to the efficiency of each module, r^j oc SiSj. 
From a combined analysis of the coincidence rates and single rates, the uncertainty 
of the extracted efficiencies was estimated to be 5%. An average coincidence rate of 
16 ± 2 Hz is observed. This agrees with the expected value of 19 ± 3 Hz, obtained 
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Figure 2: Left: Example distribution of the measured time differences between hits from neigh- 
bouring PMTs in the same storey. The sohd line is a fit to a Gaussian peak plus a flat background 
(see text). Right: Histogram of the coincidence rates (background subtracted) observed in the 
twelve detector lines. 



from a dedicated Monte Carlo simulation. The distribution of observed coincidence 
rates is shown in Figure [2] (right). The spread of the coincidence rates corresponds 
to a 10% spread in the PMT and front-end electronics efficiencies. The measured 
efficiencies have an overall normalisation uncertainty of about 15%, originating from 
the limited knowledge of the quantum efficiency and angular response of the PMT. 

3. Measurement of the atmospheric muon flux 

A muon passing through the detector can produce correlated signals in several 
storeys of the line. In the following, hits from two adjacent storeys with a local 
coincidence in each storey are considered, where the local coincidence is defined as 
a pair of hits detected by two different PMTs of the same storey within a ± 20 ns 
time window. The rate of accidental coincidence events between adjacent storeys 
is about 0.1 Hz. The data have been recorded during dedicated data taking runs, 
when bioluminescent activity was low compared to ^''K background. 

The distribution of the measured time differences between hits in adjacent storeys 
is presented in Figure [HI A clear peak is visible centred around +20 ns, demonstrat- 
ing that the majority of muons are downgoing. The width of the distribution is 
mainly due to the angular distribution of the atmospheric muons. The flat back- 
ground is due to random coincidences from and bioluminescence. The tails, 
which extend beyond ±65 ns, can be attributed to multiple muons and light scatter- 
ing. The muon event rate, R, is deflned as the integral of the peak after subtraction 
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Figure 3: Distribution of the measured time differences between hits from adjacent storeys (lower 
- upper). A local coincidence is required in both storeys. The dashed line corresponds to a Monte 
Carlo simulation based on MUPAGE (see text). 



of the flat background. On average, R is found ~ 0.06 Hz. Since the probabihty 
that a decay is seen by two PMTs separated by 14.5 m is very small, the "^^K 
contamination in R is negligible. 

The results of a Monte Carlo simulation based on MUPAGE 13j are also pre- 
sented in Figure [3l The detector response was simulated using a Geant-4 based 
model of the PMT angular acceptance and includes the effect of the detector ineffi- 
ciencies. As can be seen from Figure [3l the predicted peak is in agreement with the 
observations. 



4. Depth dependence of the atmospheric muon flux 

The rate of genuine coincidences between adjacent storeys is mainly due to the 
flux of vertical and nearly vertical muons. The muon flux includes, by definition, 
the corresponding contributions from multiple muons. The energy threshold for 
detecting downgoing muons is determined by the minimum track length to reach 
two adjacent storeys, and is about 4 GeV. 

The efficiency to detect atmospheric muons can be characterized by a single 
parameter. A, which gives the ratio between the event rate and the muon flux, 
and can be interpreted as an effective area. From the Monte Carlo simulation 
based on MUPAGE, the value for a pair of nominal storeys is A = 87^30 m^. The 
uncertainty originates from the limited knowledge of the quantum efficiency and 
angular acceptance of the PMTs, as well as uncertainties in the measurements of 
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Figure 4: The measured flux of muons as a function of depth. The results can be interpreted in 
terms of total muon intensity (left axis) as well as vertical muon intensity (right axis). The grey 
band shows the normalization uncertainty of the data. The predictions of the Monte Carlo simula- 
tions based on MUPAGE and CORSIKA are shown by dashed and dash-dotted lines, respectively. 



light absorption length in the sea water 0, 14]. A low energy cutoff of 1 GeV has 
been used for the flux integration. 

The energy spectrum and angular distribution of the muons change with depth. 
Therefore, A should be determined as a function of depth. From the Monte Carlo 
simulation, the difference between the bottom and top storey is found to be less 
than 1%. This difference has been neglected in the following. 

This analysis uses many short calibration runs taken between January 2007 and 
May 2008. During that period, in November 2007, the detector was extended from 
5 to 10 lines. The effective live times of the data samples taken with 5 and 10 
lines are 4 and 3 hours respectively. The two subsets of data have been merged 
together. The number of active storeys on each line and the efficiency of each storey 
is measured following the calibration procedure described above. From these the 
average efficiency of all the storeys at the same depth is computed via Monte Carlo 
simulation. The average efficiency for this method is about 58% with a systematic 
uncertainty of less than 10%. Because a coincidence between two storeys is required, 
a measurement of the event rate can be made at 24 different depths. The event rates 
are then corrected for the computed efficiencies and converted into flux units using 
the effective area A. The depth dependent muon flux thus obtained is shown in 
Figure HI 



The observed muon flux decreases with increasing depth, h, as expected from the 
energy loss of the muon in the water. An exponential function has been used to fit 
the data, = % x exp{{ho - h) / X) . We find % = 1.71 ±0.02 (stat) tfjl? (syst) x 
10^^ m~^ s~^ at ho = 2000 m, with a slope given by A = 540 ± 25 m. It is worth 
noting that the large uncertainty of A only affects the normalization of the measured 
dependence ($0)5 but not the slope (A). The uncertainty in A is dominated by the 
statistical uncertainties of the individual data points and is given by the fit. 

The values obtained from the Monte Carlo can be parametrized in the same 
way. This yields $mupage ^ ^ 54 ^ ^q-s ^-2 g-i 2OOO m and a slope A^^^^^^ = 
560 m (Figure HI). The data are in agreement with the simulation within a relatively 
large scale error. A similar result is obtained from a Monte Carlo simulation based 
on CORSIKA [15|]. The simulation incorporates the QGSJET model of hadronic 



interactions and the NSU model of the primary cosmic ray spectrum [16j. The 
propagation of muons in the water was simulated with MUSIC [Q. The simulation 
yields $corsika = 1.35 x 10-^ s'^ and A^°^^^^^ = 570 m. This result, also 
shown in Figure HI is compatible with the calculation from MUPAGE. 

The muon flux is often presented in terms of vertical muon intensity 18, lol. 
This usually implies a deconvolution of the measured zenith angle distribution. The 
method presented here does not require such a deconvolution procedure and provides 
the depth dependence of the muon flux directly. Hence, the results can be expressed 
in the form of vertical muon intensity. For this a constant conversion factor K, which 
translates the total flux to the corresponding vertical intensity at the same depth, is 
used. Using the Monte Carlo simulation, a conversion factor K = 0.58 sr~^ is found. 
The uncertainty of K is small compared to other normalization uncertainties. The 
variations with depth are negligible in the depth range of the ANTARES detector. 
The result of this conversion is also shown in Figure HI right vertical scale. Within 
uncertainties these results are compatible with the results of @, 0] , which refer to a 
broader depth range, partly including the depth range covered in this analysis. 



5. Summary 

A simple method for the determination of the atmospheric muon flux and its 
dependence on depth in the ANTARES neutrino telescope has been presented. 
The method is based on the measurement of photon coincidences between adja- 
cent storeys and has a low energy detection threshold. The atmospheric muon flux 
has been measured in the depth range from 2030 to 2380 m with a step of 14.5 m 
using a combined data sample of 5 and 10 line detector conflgurations. The data 
have been corrected for the presence of dead channels and unequal efficiencies of the 
PMTs, which were measured with a novel calibration technique using the natural 
radioactivity of sea water. A reasonable agreement is found between our data and 
Monte Carlo simulations. 
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